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Abstract: Nonlocal stringy versions of the Nambu-Jona-Lasinio and Gross-Neveu models 
arise in a certain limit of holographic QCD. We analyze the phase structure at finite 
density and temperature at strong coupling in terms of probe branes in the gravity dual. 
Comparison with the phase structure of the local field theory models shows qualitative 
agreement with some aspects, and disagreement with others. Finally, we explain how to 
construct the Landau potentials for these models by taking the probe branes off-shell. 
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1. Introduction 

Although most often studied as a candidate theory of quantum gravity, string theory has 
a rich history as a tool for investigating the dynamics of strongly coupled quantum field 
theories. The real-world application of most interest is QCD, an infra-red confining gauge 
theory with a small number of colors and flavors. Although the description of QCD in 
string theory remains elusive, there has been significant progress. A key development was 
the realization by Maldacena O that the near-horizon dynamics of a large number of 
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D3-branes can be expressed in two different ways, resulting in the famed duality between 
M = 4 super Yang- Mills and IIB string theory on AdS§ x S* 5 . Since that insight, systems 
of branes have been used to construct gravitational duals to field theories with flavor 0], 
with reduced supersymmetry, and that exhibit confinement /deconfinement transitions ||. 
These developments have brought us closer to capturing the dynamics of real-world QCD 
within string theory. 

One natural target for such investigations is to use string theory to learn about the 
phase diagram of QCD at finite temperature and baryon number. The QCD phase diagram 
has received much attention in recent years, and a remarkably rich and intricate structure 
has emerged, depending sensitively on the number of active flavors of quarks in the theory; 
see Q for a review and further references. Possible applications include the physics of 
heavy ion collisions and the structure of neutron stars. At the same time, the theoretical 
tools for investigating these matters are quite limited. Lattice methods are powerful for 
studying finite temperature QCD, but the notorious sign problem puts the high density 
regime out of reach. At asymptotically high density QCD becomes weakly coupled and 
tractable, but in the intermediate density regime relevant to the real world there exist no 
controlled calculational methods. 

A popular avenue of attack is to use modern variants of the Nambu- Jona-Lasinio model 
Q to study the formation of various quark condensates (important in, for example, chiral 
symmetry breaking) . The main reason why N JL models are relatively tractable is that they 
lack gluons and hence confinement. While this obviously limits their application to many 
aspects of QCD, at high density confinement does not play a dominant role {e.g. color can 
be spontaneously broken), and so NJL models are useful phenomenological models (see 
H 0] for reviews of the NJL model applied to QCD). 

More recently, an embedding of NJL-like models in QCD has been suggested in (8|. 
By taking a decompactified limit of the model of Sakai and Sugimoto M, the authors of 
Q considered a D-brane system that cleanly exhibits chiral symmetry breaking (see also 
|i"0| , |TT| ) . The brane construction of |^] provides a gravitational dual to a theory we will refer 
to as the stringy NJL model. Although the precise dynamics of this theory are unknown, 
the nomenclature follows from the appearance of chiral symmetry breaking and the lack 
of confinement. The setup involves a large number of D4-branes which intersect a stack 
of D8-branes and a stack of D8-branes, these stacks being separated in a direction along 
the D4-brane worldvolume. In the regime in which the number of 4-branes (the "color" 
branes) is parametrically larger than the number of 8-branes (the "flavor" branes), the 
strong coupling dynamics of the field theory can be easily studied using the DBI action of 
the 8-branes embedded in the near horizon background sourced by the 4-branes. 1 Unlike in 
the Sakai-Sugimoto model, chiral symmetry breaking is not explicit in the decompactified 
theory but instead emerges dynamically as the recombination of the unstable brane/anti- 
brane system, relating the separate U(Nf) factors associated with each stack of 8-branes. 

The main purpose of the present work is to map out the phase diagram of the stringy 

*A basic requirement for studying finite density QCD in string theory is to include matter in the funda- 
mental representation of the gauge group. For Nf/N c <C 1 this can be accomplished by embedding probe 
branes in a supergravity background H ; for a sampling of further work in this direction see JLa |l3[ uM . 
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NJL model by considering the theory at finite temperature and charge density. In the string 
description, this is understood in terms of exciting the gauge fields that propagate on the 
8-brane worldvolume. This topic has been studied in [l6|] 2 , but we examine this problem 
anew with an eye towards some important aspects not considered previously. In particular, 
since the NJL model lacks confinement there are finite energy massive quark states that 
are expected to play a role in the phase diagram. To incorporate this, we find a new brane 
solution that describes a condensate of smeared fundamental strings stretching from the 
8-branes to the horizon of the 4-brane geometry. Additionally, we study the possibility 
of non-homogeneous phases in which the field theory is in a configuration where regions 
of charged vacuum and uncharged vacuum coexist. It turns out that such mixed phases 
dominate the thermodynamics. We map out the resulting phase diagram and compare 
with an analogous phase diagram of a field theory NJL model. 

There are of course many limitations and caveats in attempting to apply this approach 
to real QCD. Among these, we need to be able to move away from the N c ^> 1, Nt <C N c 
regime, to incorporate bare quark masses, and to better understand the extrapolation from 
strong to weak coupling. These are all serious obstacles, but given the richness of finite 
density QCD combined with the scarcity of reliable methods for its study, we feel that 
there is more insight to be gained from developing the stringy models considered here. 

While the relation between the D-brane constructions discussed here and the local 



NJL model is quite subtle (see [11 1 for discussion), the situation is much clearer for its 



1 + 1 dimensional cousin, the Gross-Neveu (GN) model |19|. The stringy version of the 



GN-model was constructed in [10] as intersecting D4 and D6 branes. The study of this 
theory at finite charge density is especially simple. In field theory, bosonization reduces the 
problem to one involving a free scalar field. We derive the analogous result in the string 
version, which involves the dynamics of pure gauge degrees of freedom. The massless 
Goldstone excitation which relates to the phase of the chiral condensate in the field theory 
description is mapped to a pure gauge potential on the probe D6 brane whereas the massive 
fermions correspond to gauge potentials that are sourced by fundamental strings ending on 
the probe branes. As a consequence of the Chern-Simons term on the 8-brane they do not, 
however, induce charges on the non-local GN-model at the boundary of the AdS-space. 

Finally we consider the problem of constructing off-shell Landau potentials for the the 
order parameter in these models. We show that the (necessarily not unique) Landau effec- 
tive potential for the order parameter can be constructed by taking the probe-brane profile 
off-shell. In the GN model this potential can be used to reproduce the "old fashioned" 
phase diagram of the GN-model in the mean field approximation while in the D4/D8- 
brane model it can be used to study the phase structure at zero temperature. 

This paper is organized as follows. In section 2 we briefly review the .D-brane construc- 
tion of H leading to the stringy NJL model. In section 3 we consider the strongly coupled 
regime of this model, described by supergravity. We find solutions of the probe 8-branes 
with non-trivial worldvolume gauge fields. We compare the actions of these solutions in 
section 4, and so map out the phase diagram of the model. In section 5 we include a short 



2 Recent paper addressing the same question in the Sakai-Sugimoto model are jl7|, 
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but thorough review of the thermodynamics of the NJL field theory, deriving its phase 
diagram. We compare the phase diagrams of the field theory and supergravity descriptions 
in section 6. In section 7 we comment on the role of the 6-brane gauge fields in the stringy 
Gross-Neveu model. The construction and analysis of off-shell Landau potentials for the 
stringy NJL and GN models are discussed in section 8. We close with some comments in 
section 9. 



2. Weakly coupled D-brane description 

Here we briefly review the NJL D-brane construction of || (for the stringy Gross-Neveu 



model see [ 10 1 ) . In the limit of weak string coupling, the system is composed of three 
separate stacks of Z)-branes in M > . There are N c coincident L>4-branes extended in the 
{xo, X\, X2, X3, X4} directions. 3 Furthermore, there are two stacks of coincident Z)8-branes 
and coincident D8-branes, each with Nf branes. The 8-branes are extended along all 
directions but X4; in the X4 direction the stacks are separated by some distance L which we 
take to be much larger than the string length l s . We note that this configuration preserves 
no supersymmetry. 

We will study this system in the limits of low energy and weak string coupling. We 
consider a finite number Nf of flavor 8-branes and an infinite number of color branes, 
taking the 't Hooft limit of the 4-brane system, i.e. g s — > and N c — > 00 with g s N c held 
fixed and finite. In this regime, the 't Hooft coupling on the 8-branes vanishes, and so the 
U(Nf)i x U(Nf)n gauge symmetry becomes an effectively global symmetry. 

There are two (3 + l)-dimensional intersections of the 8-branes with the D4-branes. 
At low energy, the spectrum of open 4-8 strings consists only of massless chiral fermions 
localized at each intersection, left-handed for the D8-branes and right-handed for the D8- 
branes. All the fermions are in the fundamental representation of U(N C ), while the left- 
handed (right-handed) fermions are in the fundamental of the U(Nf)L (U(Nf)n) factor of 
the global flavor symmetry and are uncharged under the other factor. 

Fermions of opposite chirality interact via exchange of the (4 + l)-dimensional U(N C ) 
gluons of the D4-brane worldvolume across the gap between the stacks of 8-branes. In 
the limit of weak 't Hooft coupling (g s N c <C 1) we imagine integrating out the (4 + 1)- 
dimensional gluons to yield a theory of pure quarks where the U(N C ) symmetry is now 
understood to be global. The precise dynamics in this limit is not known and we will use 
a standard local NJL quantum field theory as a toy model. In the limit of large 't Hooft 
coupling, the dynamics are best captured by supergravity. Specifically, the 4-branes source 
a gravitational and 4-form magnetic flux background. The low-energy physics is described 
by the DBI action of the 8-branes in the near-horizon region of this background, which we 
discuss presently. 



3. Strong coupling analysis of the stringy Nambu-Jona-Lasino model 

The strong coupling description of the finite-temperature stringy NJL model with N c ^> Nf 

3 In the Sakai-Sugimoto model [B] the £4 direction is compactified. The present configuration can be 
considered a limit of that model in which the confinement scale of the dual field theory is taken to zero. 
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is encoded in the classical dynamics of Nf probe 8-branes in the near horizon geometry of 
N c non-extremal 4-branes. The worldvolume of the 4-branes lies along the Euclidean time 
r and the Xi directions, where i = 1 . . . 4. The geometry and dilaton are given by 



ds 2 



15 [f(U)(drr + (dx,y + 



+ U 2 ds% 



(3.1) 



where the non-extremality of the 4-brane background is encoded in the function 

f(U) = l-^t. (3.2) 



The horizon of the black brane is located at U = Ut- The absence of a conical 
singularity at the horizon determines the periodicity of Euclidean time r ~ r + j3, and 
hence the temperature, with 

1 
T 



ft 



47^2 



3U 



1/2 



(3.3) 



T 



and hence f(U) = 1. 



anti-D8 




U=U T 



U=U 



The limit of zero temperature corresponds to Ut 

As discussed in the previous section, chi- 

ral and anti chiral fermions are introduced D8 
by placing flavor D8 and Z)8-branes in the 
D4-brane background Jl4| , |9[ . There are two 
qualitatively distinct configurations of D8 branes. 
First, there are straight D8 and D8-branes 
separated in X4 direction, which extend all 
the way to the horizon at U = Ut- For 
the straight branes the U(Nf) x U(Nf) chi- 
ral symmetry is unbroken. Second, the D8- 
brane and anti-brane stacks can recombine 
and produce a single curved D8-brane stack 
which does not reach the horizon. For the 
curved brane the chiral symmetry is broken 
to U(N f ) §. 

A finite chemical potential for fermions can be introduced by an imaginary gauge field 
Aq = iji on the worldvolume of the 8-branes 16]. We will discuss the straight and curved 



u=oo 



Figure 1: Probe D8 branes in the DA brane 
background with attached F-strings. 



brane solutions with chemical potential in the following. One important observation is that 
it is necessary to introduce fundamental strings ending on the 8-brane to obtain non-zero 
charge density in the chirally broken phase. 

3.1 Straight brane 

The straight (anti)-D8 branes have a fixed X4 position and lie along the t, x±, X2, X3 as 
well as the U and the S 4 directions, which we use to parameterize the worldvolume. We 



4 We use units in which a' — 1. 
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choose the D8-brane stack to be localized at X4 = —L/2 and the anti-Z)8-brane stack at 
X4 = L/2. The worldvolume of the branes reaches the horizon at U = Ut- To include 
a chemical potential in the dual description, we turn on the 8-brane gauge field Aq and 
assume all other components vanish. The Dirac-Born-Infeld action for the D8-brane as 
well as the D8-branes is then given by 



00 



S s = kR^ dUU^ J 1 + (A' ) 2 , (3.4) 
Ju T 

where primes denotes derivatives with respect to U. The constant coefficient k is given by 

T 

K = —N f {3V 3 n 4 , (3.5) 
9s 

where Tg is the 8-brane tension and V3 and are the volumes of the three noncompact 
directions {£1,2,3} and the unit S , respectively. 

Varying the action (|3.4| ), we find an integral of motion (constant in U) 

iQs = ^r = ^ UlA '° ■ (3-6) 



5A' 







1 + (A' ) 



2 



We will call Q s the "conserved charge", although it is more properly the electric flux on 
the brane. We can integrate (^^) to solve for the gauge field. This relates the conserved 
charge and chemical potential /j, = IAq{U = 00) via 

t i = Q s [°°dU - 1 , (3.7) 

Ju T y/K 2 R 3 U 5 + Qf 

where we have set Aq(Ut) = to avoid a singularity at the horizon. 

The combined action of the branes and anti-branes in terms of the conserved charge 
is given by 

roo tt5 

S s = 2kR 3 / 2 / dU — . (3.8) 

Note that both the brane separation L as well as the non-extremality function f(U) are 
irrelevant for the straight brane action. For the straight brane as well as the curved brane, 
the Chern-Simons term J AAdAAF§ on the D8-brane worldvolume can be neglected, since 
the only nontrivial gauge field component is Aq. 

3.2 Curved brane 

The brane configuration corresponding to broken chiral symmetry is a curved Z)8-brane. 
We treat this similarly to the straight brane configuration except that we use the coordinate 
X4 instead of U to parameterize the worldvolume. The embedding is determined by the 
profile function U(x±) with boundary conditions such that U{x± — * ±L/2) — * 00. The 
curved brane is symmetric about X4 = 0, where its position f/(a?4 = 0) = Uq is closest to 
the horizon. As we shall see it is necessary to introduce fundamental string ending on the 
Z)8-brane to achieve a nonzero charge density. 
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The Dirac-Born-Infeld action for the D8-brane in this parametrization is 



S c = « dx, U i Jf(U)+(^ \u 2 + Al) , (3.9) 

where the dot denotes a derivative with respect to x±. There are integrals of motion 
independent of X4 



P 4 t Si) i c - U ' f W (3 . 10 ) 

Once again we refer to these integrals of motion as conserved charges. In order to obtain 
the same chemical potential on both branches of the D8 brane the gauge field Aq must be 
symmetric about X4 = 0. This implies that, for a smooth gauge field, Aq = at X4 = 0. It 
follows from ( |3.10| ) that Q c = and hence Aq has to be constant. Since the Chern-Simons 
term is not relevant for the stringy NJL model (in contrast to the stringy GN model, as 
we will see later) the constant gauge field can be set to zero and hence it seems that the 
chemical potential has a trivial effect on the dynamics. 

A way to overcome this is to loosen the condition that the gauge field Aq and the 
brane embedding U(x4) must be smooth at X4 = 0. This can be achieved by introducing 
fundamental strings located at X4 = stretching from the horizon to the worldvolume of the 
Z)8-brane resulting in the configuration shown in Figure ||. The action for Ni fundamental 
strings is (recall we set a' = 1) 

= ^ Jdx A [U(x 4 ) - U T ]S(x 4 ) + % *^-jdx A 5(x 4 )A . (3.11) 

The variation of the combined action of fundamental strings and D8-branes leads to the 
jump conditions for the Q c and U at X4 = 

AQ C = 2Q C V ' 



rt-rt^fMf*)*™, (3.12) 



2tt 

/ 

p 4 V R J 27r 

where we used the fact that for the symmetric L>8-brane the gauge field derivative Aq and 
U change sign at X4 = 0. We choose the charge on the right side of the D8-brane (where 
X4 > 0) to be positive. Using ( p.!2| ), together with ( p. 10 ), one can show that U = —iAq at 



X4 = 0. Substituting this back into the expression for P4 given in ( 3.10 ) one finds 



Pa = KU$Vf(U ) ■ (3.13) 
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The expression of the conserved charges ( 3,10| ) can be solved for the derivatives Aq and U{ 



o 



o _ n%_ (<L\ 5 .La. + C^V « _ g» 



taJiVWd w V w /(c) 

_M„ _ * . (3.14, 

Integration of the first equation in (|3.14| ) relates the conserved charges to the width sepa- 
ration L in the X4 direction of the two asymptotic legs of the D8-brane. 



L 



00 dU 
u ~U 



, r°° r iTj 

= 2C^v7(^o) / 5 , ■ (3-15) 

Integration of the second equation in ( |3,14| ) relates the conserved charges and the chemical 
potential 

rco ■ X 

^ = -iA {U )+ dU " 



u U 

iA (U Q ) + ^ f°° 3 ; dU ■ ( 3 - 16 ) 



These two equations relate the parameters governing the curved brane: L, fj,, Q c , Uq and T. 
For example at fixed temperature, ( p. 15 ) and ( 3.16 ) can be used to find the turning point 
Uq and the charge Q c , given the separation L and the chemical potential fi and vice versa. 



4. Phase diagrams at strong coupling 

The solutions discussed in the last section provide strong coupling descriptions of possible 
phases for the stringy NJL model. To determine which phase dominates for specific exter- 
nal conditions, we must compare their free energies (or other appropriate thermodynamic 
potentials). In the language of the previous section, this entails comparison of on-shell 
actions (or appropriate Legendre transforms). 

When including the effects of non-zero gauge field on the probe 8-branes, one can 
perform the comparisons at fixed charge density or fixed chemical potential. We will 
analyze the two cases in turn. 

4.1 Fixed chemical potential 

The thermodynamic potential which governs the phase structure at fixed chemical 
potential is £l(T,n) = E — TS — fiQ, the grand free energy. In the semiclassical probe 
approximation the grand free energy of the system at fixed chemical potential is given by 
the Dirac-Born-Infeld action of the branes and the Nambu-Goto action of the fundamental 
string, where the boundary value of the imaginary gauge field Aq(U = 00) = i[i determines 
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the chemical potential. The actions are formally divergent so in order to evaluate them we 
introduce an ultraviolet cutoff at U = A. All actions have the same divergence in A and 
so the difference of the grand free energies is independent of the cutoff as A — ► oo. 

In the following we choose units of length s_c(qj 
and action such that R = 1 and k = 1. The 
grand free energy for the unbroken phase is 
given by the action of the straight branes 
(13. 8f) and reads 



A z 5 
dz — ■ 



-A 7 / 2 
7 



(4.1) 

The charge density Q s on the straight branes 
is determined uniquely by the chemical po- 




tential through (3.7). 

The grand free energy for the broken 
phase is given by the sum of the action of the curved brane ( 
mental strings ( |3.11 



Figure 2: Action of the curved brane as a func- 
tion of the charge Q at T — .3 and L = .5 for 
various values of the chemical potential /i 



and the smeared funda- 



S c (v,Qc) = 2 



dz 



U 



U§f(Vo) 



+ Q c (Uo - U T ) + iA (U )Qc 



2 A 7/2 



(4.2) 



Unlike the straight brane scenario, the charge density is not determined uniquely in terms 
of the boundary data (the chemical potential /i and the asymptotic brane separation L). 
Rather, even after imposing the constraints ( |3.15 ) and (3.16) to eliminate the location Uq 
of the tip of the curved brane and the value of the gauge field at the tip Aq(Uq), the charge 
is still free to be varied. The charge on the curved brane is then that which minimizes 
the ( f4.2|) for a fixed chemical potential. The behavior of the action displayed in figure || 
is generic; for small chemical potential the Q c = curved brane is dominant whereas for 
larger chemical potential the curved branes with Q c 7^ are dominant. 

To obtain the phase diagram of the system, we rely on numerical comparison of the 
actions discussed above. We find that at sufficiently low temperature, the curved brane 
without charge density is thermodynamically preferred and so chiral symmetry is broken. 
Raising the temperature and/or potential, one encounters a transition to the straight brane 
(i.e. chiral symmetry is restored), whose charge is determined uniquely by the chemical 
potential. Although the Q c 7^ phase of the curved brane does eventually become preferred 
over the Q c = phase, it appears that it never is preferred over the straight brane phase. 
We show the phase diagram for the system in figure ||. 

The phase boundary is given by the curve in //, T space where £ c (/i, Q c = 0) = S s (n). 
The phase diagram displayed in figure ||| shows this phase boundary, the critical value of 
the chemical potential fiphase as a function of the temperature T. At a critical temperature 
of T = .336 the curved brane ceases to exist and the system is in the unbroken phase. The 
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plot is presented for a specific value of the separation L. However as we shall argue in the 
end of this section, a scaling argument shows that the phase diagram is universal. 



4.2 Fixed charge 



In the previous section we worked at fixed 
chemical potential, where the thermodynamic 
potential which governs the phase structure 
is f2(T, /i) = E — TS — ^jlQ. If instead we work 
at fixed charge the new thermodynamic po- 
tential is the Helmholtz free energy, given by 
a Legendre transform F(T, Q) = £1 + fiQ. 
On the gravity side this is accomplished by 
adding fiQ = —iAo(U = oo)Q to the action. 
Note that for a given fixed charge, the gauge 
field Aq(U = oo ) and hence the chemical po- 
tential [i is not the same for the straight and 
curved brane, but is determined in terms of 
the charge by the constraints (3.7) and ( |3.16|) 
respectively. 

The modified action for the straight brane is given by 

S,(Q) = S s + ix s Q = 2( ( K dUy/U^ + Q 

v Ju T 



Figure 3: Phase diagram at fixed chemical po- 
tential. Below the plotted curve, chiral symme- 
try is broken. Above the curve, the symmetry 
is restored. 



7 



-A 7 / 2 



The modified action for the curved brane with fundamental strings attached 



S C {Q) = S c + ijl c Q = 2 



dz 



z 5 + Q 2 



U 



+ Q 2 



+ Q(U - U T ) 



-A 7 / 2 
7 



(4.3) 



(4.4) 



As before the phase diagram can be obtained by considering the sign of AS = S s — S c . 
The phase boundary in the Q,T plane for separation L = 0.5 is plotted in figure ||. It is 
an interesting feature of the phase diagram that the critical value of the charge increases 
as the temperature is increases. 

The numerical analysis of the phase structure in the last two sections was done at a 
fixed value L = 0.5 of the separation. The phase structure is however universal in the 
following sense. If we scale L — ► La along with 



U^U/a 2 , Q -» Q/a 5 , n^fi/a 2 



(4.5) 



the constraints (p~7|), ( |3.15 ) and ( |3.16| ) are invariant. Under these scalings the actions 
Q4.3| ) and ( |4.4| ) both scale as S — > S/a 1 . Since the phase structure is only sensitive to 
the relative sign of the difference of actions the phase diagram remains the same under 
the scaling. Although the solutions depend in general on three parameters L,T, \x at fixed 
chemical potential and L, T, Q at fixed charge, the scaling symmetry can be used to fix one 
of them, which we have done. 
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4.3 Mixed phase 

In the case of fixed charges we considered 
two phases: the straight branes and curved 
brane with constant charge density. As dis- 
cussed earlier the curved brane must have 
fundamental strings ending on it, in order to 
support a nonzero charge density. 

There is however the possibility to have 
a inhomogeneous mixed phase with regions 

of broken chiral symmetry and regions where 

Figure 4: Phase diagram for fixed charge with 
the symmetry is restored. The charge is car- £ — q 5 

ried by the straight brane regions and the 

curved brane regions do not have any fundamental strings ending on them. If the total 
charge density is Q and the straight brane takes up yV of the volume V whereas the curved 
brane with no charge takes up (1 — y)V of the volume, the charge density on the straight 
brane has to be Q/y. If the contribution of the interface energy is negligible (which seems 
sensible in the large volume limit since the bulk contribution scales as the volume whereas 
the interface scales as the area) the free energy in the mixed phase is given by 

*"«■'>- <*£*f+® + (i -<^7j4pgH A7/2 ) ■ 

(4.6) 

The free energy of the mixed phase should be compared to the free energy of the curved 
and straight brane at the same total charge density. Numerically we have investigated this 
question and we checked that for all our choices of charges and temperatures the mixed 
phase has smaller free energy than the curved brane with the same charge. 

Hence for parameters where the previous analysis showed that the chiral symmetry is 
broken, the preferred phase of the system is mixed. Incidentally for some range of charges 
this is also true for the phase where the straight brane dominates since the free energy of 
the straight brane is given by setting y = 1 for the mixed phase. For small charges even 
then the mixed phase dominates. If the charge is increased the straight brane will have 
smaller free energy and dominate. A representative plot of the free energy as a function of 
y for various values of the charges is given in figure [|. Note that at large enough charge 
the minimum free energy is at y = 1, indicating a homogeneous phase of restored chiral 
symmetry. 

5. NJL model at finite temperature and chemical potential 

In the previous section we studied the finite temperature/chemical potential phase diagram 
of the stringy NJL model at strong coupling. Now we want to compare against an analo- 
gous phase diagram in ordinary field theory. The field theory in question will be the usual 
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NJL model at large N. Before turning to the calculations we pause to motivate our study 
of the NJL model in the present context. As noted in the introduction, the eventual goal is 
to learn something about finite density QCD. The phase structure of high density QCD is 
for the most part governed by the condensation of various fermion bilinears. For instance, 
at sufficiently high density it is known that a diquark condensate forms, spontaneously 
breaking the gauge symmetry ( "color superconductivity" ) . The formation of such conden- 
sates can be studied in the context of NJL models with four-fermion interactions. Indeed, 
at asymptotically high density such a treatment is under good control due to asymptotic 
freedom. At the smaller densities relevant to the real world, controlled computations are 
harder to come by, but NJL models still serve as valuable phenomenological guides, al- 
though one should be aware of ambiguities due to cutoff dependence and the neglect of 
higher order fermion interactions. 

The strong coupling stringy NJL model sjnued 
studied here is an example of a specific NJL- 
type model, whose virtue is that it can be 
easily studied using holography. Ideally, we 
would like to be able to do a corresponding 
field theory computation of the phase struc- 
ture of the same model. Instead, we just 
compare against the simplest large N NJL 
model, which should therefore be thought of 
as a toy model of the real problem. Also, we 
will only study the formation of ipip conden- 
sates, even though for real QCD condensates 
of the form ipip are relevant for color super- 
conductivity. Our reason is simply that in 

our supergravity description we only saw the appearance of ipip condensates. A supergrav- 
ity model of color superconductivity would of course be interesting to study. 

We now turn to a review of the finite temperature/density phase diagram of the large 
N NJL model. This model can be studied using standard large N techniques. Some 
relevant references are [20|, 21, |22| , 0] . 

The general NJL model is a theory of Dirac fermions if) za , with i = 1 . . . Nf a flavor 
index, and a = 1 . . . N c a color index, that exhibits spontaneous breaking of a continuous 
chiral flavor symmetry. A particular example is 



Q=2.0 



Q=0.6 



Figure 5: Free energy versus fraction y of 
system in chirally symmetric phase. Plots for 
L = 0.5, T = 0.2 and three distinct charges 
Q = 0.2, Q = 0.6, Q = 2.0. 



£ = m)i> + G [(%bT A iP) 2 + (iPi-f 5 T A %by 



(5.1) 



where T A are generators of U{Nj). Roughly speaking, we think of the four-fermi interac- 
tions as being induced by single gluon exchange between quarks. Of course, in real QCD 
there is no justification for omitting multiple gluon exchanges, which correspond to higher 
order fermion interactions. A related point is that the four-fermi terms in ( |5.1| ) render 
this theory nonrenormalizable, and so a UV cutoff A needs to be introduced to define the 
theory. While the precise quantitative predictions derived from this model are strongly 
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cutoff dependent, the NJL model seems to accurately model the qualitative effects of chiral 
symmetry breaking in QCD. 

For sufficiently large G this theory develops a nonzero fermion bilinear condensate. 
The mass term m represents explicit breaking of chiral symmetry, and yields a preferred 
orientation of the condensate, namely (ipip) 7^ 0, leaving an unbroken diagonal U (Nf) flavor 
symmetry. We then think of taking m — > to remove the explicit breaking; the motivation 
for this in the present context is that the corresponding mass term is absent in the brane 
picture. 

At large N = NfN c the standard approach is to introduce an auxiliary field for the 
condensing bilinear, 

£ = 2VG^ - (f> 2 + G [(W A ipf + $iKT A iP) 2 ] . (5.2) 

Here the notation T' A indicates that the U(l) generator is omitted. Upon integrating out 
the auxiliary field <j> we reproduce (5.1). 

If we instead integrate out the fermions we induce an effective potential for eft. The 
expectation value of (f) yields the fermion condensate via = (4>) ■ The leading large 
N contribution to V((f>) comes from 1-loop diagrams with a fermion in the loop and an 
arbitrary number of external (ft insertions. From the structure of ( |5.2D we see that at this 
order we only need the action 



C = VX^+ 2VG4>)i) - <p 2 . (5.3) 
The effects of the omitted terms are sub leading in 1/N. 
5.1 Zero temperature and chemical potential 



Integrating out the fermion in fl5.3| ) in the presence of constant <p yields the effective po- 
tential 

V((j)) =cp 2 + iTi lnfof + M) , (5.4) 
where the effective fermion mass M is 

M = -2\fG<j) . (5.5) 

The "gap equation" is obtained by minimizing with respect to (j). This gives 

t-^Ta—o 7775 • (5.6) 

(2vr) 4 p 2 — M 2 + it y ' 

Wick rotating and evaluating the momentum integral with a hard cutoff, \p\ < A, we obtain 



1 = A 



fA 2 Y\ ( A 2 

{w) ln { 1 + w 



(5.7) 



where we defined the dimensionless coupling 



A = . (5.8) 
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To determine the condition for chiral symmetry breaking we note that f(x) = 1 — 
~ ln(l + x) is a monotonically increasing function, obeying /(0) = and /(oo) = 1. Thus 
the gap equation (5/7) has no solution for A < 1. Since in going from fl5.6Q to ( |5.7D we 



divided by M, this implies that for A < 1 we have M = 0, so that chiral symmetry is 
unbroken at weak coupling. For A > 1 the gap equation admits a solution with M/0, 
and it is easy to check that this minimizes V{4>)- Hence chiral symmetry is spontaneously 
broken for A > 1. 

5.2 Finite temperature and chemical potential 

We first recall a few basics. The grand partition sum is 

Zn = Y,e-^ Q) , (5.9) 
and its logarithm gives the grand free energy 

0(T,/i) = -ilnZn = E-fiQ-TS . (5.10) 

At fixed T and \i the preferred phase of the system is the one that minimizes O. If we 
instead work at fixed charge Q, then we should minimize the Helmholtz free energy , 

F(T,Q) = E -TS = Q + \iQ . (5.11) 

Returning now to the NJL model, we want to find the gap equation at finite tempera- 
ture and chemical potential. It is straightforward to generalize the previous derivation (see 
|20| , [j], |2^|), but we instead give a simple physical argument. If we go back to ( |5,6D and 
perform the p° integral we get 

M = iN < N ° G lw?w p ' <5 - 12) 



where E p = -s/p 2 + M 2 . There is an intuitive way to think about ( 5.12j ): it is the statement 



that the assumed effective mass M is indeed the mass obtained by the fermion interacting 
with the Dirac sea. This picture makes it clear that to go to finite chemical potential and 
temperature we instead write 



M = ANfN c G J ^3 ^ (l " n((3, /x) - n(/3, /.)) , (5.13) 

with 

The thermal gap equation (|5.13| ) now says that the fermion interacts with the Dirac sea 
and the Fermi-Dirac distribution of fermions and anti-fermions present. This is the same 
result obtained from a more systematic derivation. 
In the same spirit, the grand free energy is 

m /*) = ^ - ™ F N C J [ Ep + i In (l + e-^-")) + 1 In (l + e~^^ 

(5.15) 
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Indeed, extremizing f2 with respect to M yields ( 5.13| ) 



To map out the phase diagram we need to minimize f2 at fixed (/?,//). For future 
reference the charge density is 

Q((3, (i) = 2N f N c J ^3 (n(/3, y) - n((3, M )) . (5.16) 

5.3 Phase diagram 

In this section we analyze the phase diagram at fixed temperature and chemical potential. 
At low temperatures we find a first order phase transition that restores chiral symmetry 
for sufficiently large chemical potential. Increasing the temperature, we eventually find a 
tricritical point beyond which the first order transition becomes second order. 

We take G sufficiently large such that at zero temperature and chemical potential 
there is chiral symmetry breaking (this corresponds to taking A > ^). There is then a 
curve fJ,ait(T) above which chiral symmetry is restored. The curve corresponds to the 
locus of points for which the grand free energies of the broken and unbroken phases are 
equal. At low temperature, if we raise the chemical potential above n c rit(T) we continue 
to find solutions to the gap equation in the broken phase. These solutions are metastable, 
since the unbroken phase is thermodynamically preferred for fi > ^ cr i t (T). Continuing to 
increase the chemical potential, we eventually find that there no longer exist solutions to 
the gap equation in the broken phase, and so the metastable region terminates. In the 
following, we will compute both the actual phase boundary as well as the metstable region. 

Introducing dimensionless variables via p = Ap , // = A/t , (3 = A -1 /? , M = AM, where 
A is the UV cutoff on the three-momentum, the gap equation becomes 

1 = 4A f dp p2 (l - n0, p) - n0, /})) . (5.17) 
Jo Jp 2 + M 2 J 



Here we are using the same symbol for A as in (5.12), although the interpretation is slightly 
different since A is now a cutoff on 3-momentum rather than 4-momentum. We can also 
define a rescaled charge density 

Q = J^§J^ = J/P P 2 A) " n0, A)) • (5.18) 
For sufficiently low temperature and chemical potential there exist solutions to the gap 



equation (5.17) with M > 0, indicating the existence of (stable or metastable) symmetry 
breaking vacua. Solutions cease to exist for sufficiently large temperature and chemical 
potential, demonstrating symmetry restoration. We will first find the boundary between 
these two regions. 

For definiteness we choose A = 1. It is useful to first consider the special cases of 
vanishing temperature or chemical potential. 

At vanishing chemical potential (and hence charge) it is easy to check that the gap 
equation ceases to have a solution for T > T c ~ .59. Further, as the temperature approaches 
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T c from below, the equilibrium value of M moves smoothly to the origin. Thus the phase 
transition is second order at finite temperature and vanishing density. 

Next consider vanishing temperature and positive chemical potential. In this case the 
distribution functions are 



n0, fx) = B (p F - p) , n0, jX) = , 

with 

p F = \J 'jp - M 2 8(/t - M) . 
Hence the gap equation and charge density are 



1 = 4 



P 



PF ^ p2 + M 2 



(5.19) 
(5.20) 

(5.21) 



Consider raising the value of fi from 0. At fx = we have an uncharged vacuum at M ~ 1.1. 
Clearly, this remains a solution until we reach fx ~ 1.1. As we take fx > 1.1 we find that 
there is no solution to the gap equation, which means that we are driven to the symmetric 
vacuum at M = 0. However, well before (at fx ~ .83) we reach this value of fx it is 
straightforward to check that the symmetric vacuum has the lower free energy. Therefore, 
there is a first order phase transition at fx ~ .83. Note that the transition is between an 
uncharged vacuum at M ~ 1.1 to a charged vacuum at M = 0. 

Finding the complete phase diagram requires a more involved numerical analysis. We 
proceed by sampling a number of discrete choices for (T, fi, M) and checking whether 
the gap equation admits a solution or not. We then locate the boundary between the 
regions with and without solutions. Given a solution to the gap equation with M > 0, 
we need to evaluate its grand free energy to see if it or the unbroken phase at M = is 
thermodynamically preferred. The phase boundary corresponds to points where the two 
free energies are equal. We show the phase diagram in the T — fi plane in figure |(| 

The upper curve corresponds to the bound- 
ary of the region where solutions to the gap 
equation exist. The lower curve is the true phase 
boundary, where the free energies are equal. In 
the temperature regime where the two curves 
are distinct there is a first order phase transi- 
tion, with points between the two curves rep- 
resenting metastable phases. This ends at a 
tricritical point when the curves meet, and for 
higher temperatures the phase transition is sec- 
ond order. 

By using ( 5.18; ) we can translate from chem- 
ical potential to charge density. The phase di- 
agram in the T — Q plane is shown in figure 0. 




0.2 0.3 0.4 

temperature 



Figure 6: Phase diagram in the T — fi plane 
There is a tricritical point where the first order and second order phase boundaries meet. 
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temperature 

Figure 7: Phase diagram in the T — Q plane. 

At low temperature we have a first order phase transition, and a corresponding coexistence 
between the broken and unbroken vacua. At fixed chemical potential these two vacua have 
different charge densities, as indicated in the figure (the smaller of the charge densities 
corresponds to the broken phase). 

For temperatures below the tricritical point, instead of fixing the chemical potential 
were we to fix the charge density then the system will be in a mixed phase, with distinct 
regions of space being filled by broken and unbroken phase. As we increase the charge 
density we convert more and more of space into the unbroken phase. 

At low temperature, the resulting picture gives qualitative agreement with QCD if we 
think of the broken phase as vacuum and the unbroken phase as nuclear matter. Thus 
the NJL model can be a useful guide to studying cold dense matter in QCD. At higher 
temperatures, though, the broken phase of the NJL model starts to acquire a nonzero 
charge density, which has no QCD analog. This is of course due to the lack of confinement, 
and demonstrates that the NJL model is a poor guide to high temperature QCD. 

6. Comparison of supergravity versus field theory phase diagrams 

We can now make a qualitative comparison between the supergravity and field theory phase 
diagrams. It turns out that at low temperature the phase structures are in reasonable 
agreement, while they differ at high temperature. 

It is easiest to make the comparison at fixed chemical potential, which amounts to 
comparing figures ^ and ||. Consider working at fixed temperature and raising the chemical 
potential from zero. At low chemical potential we are in the broken phase. As we raise 
the chemical potential we eventually enter a range in which the broken phase becomes 
metastable. On the supergravity side this metastable phase is a curved brane with no 
strings attached, hence carrying no charge. This is in agreement with what we find in 
the field theory; there the metastable phase corresponds to fi < M, which again implies 
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vanishing charge. Now raise the chemical potential further. On the supergravity side we 
find that the metastable curved brane persists but starts to carry a nonzero charge due 
to fundamental strings. This aspect differs from the field theory, where for sufficiently 
large chemical potential the metastable phase ceases to exist (there is no solution to the 
gap equation). On both sides, the true phase transition is between an uncharged broken 
phase and a charged unbroken phase (the straight brane in supergravity and the M = 
vacuum in field theory). Furthermore, the order parameter corresponding to the fermion 
mass jumps at the transition, and so the phase transition is first order in both descriptions. 
Altogether, we find that there is reasonable qualitative agreement between the two sides, 
including the existence of metastable phases, with the differences becoming apparent at 
the largest chemical potentials. 

At higher temperatures the situation is qualitatively different. The difference is that in 
field theory there is a tricritical point at which the phase transition switches from being first 
order to second order. For temperatures above the tricritical point there is no metastable 
broken phase. On the supergravity side the transition is always first order. A second order 
transition would correspond to the branes touching the horizon, but it turns out that this 
is never a solution of the equations of motion. 



7. Stringy Gross-Neveu model 



The 1+1 dimensional analogue of the NJL model is the Gross-Neveu model [19|. In this 
section we make some observations about the finite chemical potential thermodynamics 
of this theory, in the field theory and brane descriptions. As we review, the analysis is 
essentially trivial since the chemical potential just couples to a free boson. Nevertheless, it 
is instructive to compare how this comes about in the two descriptions. 

7.1 Chiral Gross-Neveu model in field theory 

The chiral Gross-Neveu model is a 1+1 dimensional field theory with four-fermi interactions 
invariant under a continuous chiral symmetry, 

^ = ^ + ^[W) 2 -(?7 5 ^) 2 ] • (7.1) 

Here the fermions are N component objects, with the indices suppressed and A is the 't 
Hooft coupling. This theory has conserved vector and axial currents. 

The interaction has the same structure as that induced by integrating out gluons. 
Indeed, after using Fierz identities we can also write the interaction term as 

where the singlet and non-singlet currents are 

f = , f A = ^T A ^ . (7.3) 
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This makes it clear that the theory separates into a U(l) piece and an SU(N) piece |23| . 
The £7(l)-factor whose conserved charge is the fermion number ("baryon number") can be 
bosonized in terms of a canonically normalized free compact boson of radius V Nir 




f = -\-r^ v dA, (7.4) 



so that fermion number is just the winding number of (j) [23, 24 1 Adding the interaction 
just changes the coefficient of the kinetic term of <f). The nontrivial part of the theory is 
an SU(N) WZW model with added term j A j A . Or, equivalently, it is the massless SU(N) 



Thirring model 23] 



The bosonization ( |7,4| ) shows that the free scalar parameterizes the rotation of the 
condensate. We do not actually have symmetry breaking since massless scalars in two 
dimensions cannot sustain expectation values due to IR fluctuations. But for large N 
these fluctuations are suppressed, so we can still think of the symmetry breaking in an 
approximate sense. More precisely, the relevant two point correlator dies off like [rrl -1 /^ 



& 

Adding a chemical potential for baryon number only affects the U (1) part of the theory, 
and so is extremely simple to analyze. In terms of the free boson, a nonzero baryon density 
is carried by a linearly varying profile, (f> = — \J~^[ Pbx 1 . Balancing the energy density 
e = ^p\ against the chemical potential term pps, we find that the baryon number is 
proportional to the chemical potential, pb = J^/U- In the fermion language this just 
corresponds to filling up the Fermi sea. 

The finite temperature thermodynamics is nontrivial since the SU(N) part of the 
theory participates. But again, the effect of a chemical potential is simple since it only 
involves the U(l) part. The complete effect of the chemical potential is captured just by 
turning on a linearly varying background profile for the free boson. 

7.2 Brane version of Gross-Neveu model 

We now consider the non-local generalization of the Gross-Neveu model realized as inter- 
secting DA and D6-branes |jio| ]. In the strong coupling limit this system has an effective 
description as a D6 — D6 probe wrapping an 5 4 in the near horizon geometry of N DA 
branes given by fl3.1|) . The setup is analogous to that giving rise to the stringy NJL 
model, the difference being that the dual field theory is now (1 + I)-dimensional rather 
than (3 + l)-dimensional, since the D8-branes are replaced by L>6-branes. 

In the field theory analysis we saw that the theory separates into a U(l) piece and 
an SU(N) piece. The U(l) piece is described by a free boson, which carries the baryon 
number of the theory via its winding number. We would now like to see how this feature 
emerge in the strong coupling description. 

There are 2+1 unwrapped dimensions of the D6 — D6 probe. The main ingredient for 
understanding the above question is the presence of a Chern-Simons term for the world- 
volume gauge field, 

S C s = ^ [ AAdA. (7.5) 



Air 
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where the integration is over the (1 + l)-dimensional intersection as well as either the U or 
X4 directions. To derive this, we start from the RR-coupling on the D6-brane worldvolume, 

5 6 = ^(2ira') 2 T 6 J A A dA A G 4 • (7.6) 

In the presence of N D4-branes we have /54G4 = 2k1 T^N. Integrating over the S and 
using T 4 T 6 = [(2vr)(2Kf )(27ra / ) 2 ]- 1 we find the result flFigp . 

Being the lowest derivative term for the gauge field on the Z)6-branes, the Chern- 
Simons term dominates at long distance. In the usual spirit of holography we can relate 
the asymptotic behavior of the bulk gauge field to the currents and external gauge fields 
in the boundary field theory. To make this precise we first need to supplement the action 



Si 



TUp with a boundary term [^. 27] 

— / y/ggTA^ . (7.7) 

OTT Jqm 



bndy " 

Recall that the brane worldvolume boundary has two distinct components; denote by A^ 1 ' 2 ^ 
the value of the gauge field on the respective boundary component. With the above choice 
of boundary term the on-shell variation of the action is 



5S 



J d 2 x{A^5A W + aPsA®) , (7.8) 



where x^ are lightcone coordinates along the boundary. From this expression we can read 
off the currents 

Note that the left and right moving currents are each localized on one of the boundary 
components. This corresponds to the fact that the fermions localized at each intersection 
are chiral (and have opposite relative chirality). The preceding formulae hold whether we 
are in the unbroken phase (straight branes) or broken phase (curved brane). 

Next, we discuss the interpretation of gauge transformations on the D6-branes. In the 
unbroken phase we can choose the gauge Ajj = on the two straight branes, which still 
leaves us with the freedom to perform U independent gauge transformation on the brane 
and anti-brane, corresponding to an unbroken U{V)v X U(1)a chiral symmetry. In the 
broken phase we can choose A^ = along the curved brane. We then have a single function's 
worth of x 4 independent gauge transformations, which we identify with an unbroken U (l)y 
symmetry. 

The spontaneously broken U(1) A transformations can be defined as the gauge trans- 
formations 

SA = dA A , A A = f(x 4 )X A (x + , x~) , (7.10) 

where /(±L/2) = ±1. These symmetries are spontaneously broken in the sense that they 
fail to preserve A± = 0. 

We expect U(1)a to act as a shift of the associated Goldstone boson <f> (the "pion"). 

Indeed, if we define <f> as e^J~^^ = A ' d \ where the integration contour connects the two 
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boundary components at a given x^, we see that U(\)a acts as 5(f) = y Furthermore, 
is invariant under U(l)y, as it should be. This identification of the Goldstone boson 
agrees with that in |l4j ]. 

The pion <p is identified with the bosonizing field appearing in ( |7,4| ). To see this, note 
that upon acting with U(1)a we have the currents 





U = -\l— d+t, j- = yl—d-0, (7.11) 

which agree with (|7ip. 

As in the field theory, baryon number is carried by the winding number of the pion. 5 In 
the brane setup we can turn on baryon number by performing a large U(1)a transformation. 
The choice Xa = \J~^Pb induces the baryon number density ps- The energy cost for doing 
this is entirely governed by the boundary term ( |7.71 ), since this is the only nonvanishing 
contribution to the action in the presence of a flat connection. The stress tensor following 
from this boundary term is that of a free boson, in agreement with the field theory. 

It is now easy to check that the effect of a chemical potential in the brane setup will 
agree with that in the field theory. A chemical potential p is mapped to the boundary 
conditions A_ = A^' = —^fi. The easiest way to satisfy this is to just take the constant 
gauge potential = A + = — Then from ( |7.9| ) we find that ps = j+ + j- = 
which agrees with what we found in the field theory. 

To summarize the results of this subsection, we find that the results in field theory have 
a simple analogue in the brane setup. In field theory, bosonization reduced the problem 
of finite baryon density to a problem involving a free scalar field. On the brane side, 
everything reduces to pure gauge configurations. Because the branes have a boundary, 
these pure gauge configurations describe physical degrees of freedom. In particular, we saw 
that these degrees of freedom precisely match those of the bosonizing field. 

8. Landau potentials 

In our supergravity analysis of the phase structure of the stringy NJL and GN models we 
have looked for solutions in which there is a smeared distribution of fundamental strings 
attached to the probe branes. A point that we have not emphasized so far is that for 
generic values of the fundamental string charge the homogenous solution for the curved 
brane profile, and therefore the value of the order parameter for the chiral symmetry 
breaking, is not unique. In field theory near the critical point this type of behavior is 
typical in the Landau theory of phase transitions. In this section we address the question 
of how a Landau potential can be obtained in the supergravity description. At first sight 
one might conclude that this information is not available on the supergravity side since 
the correspondence between supergravity and field theory is on-shell in the bulk. 6 In 

5 This is a simpler version of the statement that in the Skyrme model of 3 + 1 dimensional QCD a 
baryon can be identified with a winding configuration of the Goldstone bosons. In that case, there is a dual 
description in terms of an instanton on the brane [ pi] , p9[ . 

6 See, however, |3(| for an example where agreement between a certain function evaluated in the CFT 
and a family of AdS - black holes is found off-shell. 
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general, it is easier to take open strings off-shell than closed strings, since the latter include 
gravity. In the present case, we can certainly take the profile of the probe brane off-shell, 
while continuing to solve the closed string equations of motion, and the free energy of this 
configuration should have a field theory interpretation. 

The strategy which we will follow is to take the brane profile off-shell by giving up 
the jump condition for the profile function U{r) at the tip of the brane. As we will see 
this leads to an effective potential for the order parameter Uq which interpolates between 
different solutions. This gives a convenient picture of the phase structure, and furthermore 
allows one to assess the (in)stability of the various extrema. In the following we will just 
be considering the zero temperature case for simplicity, but the method is straightforward 
to extend to finite temperature. 

8.1 Gross-Neveu model 

Let us start with the stringy version of the Gross-Neveu model at finite density. The 
Landau potential in this theory is somewhat artificial since there is no phase transition 
here as a function of baryon density as explained in p4| , [31]]. On the other hand, we can 
think of studying the phase structure in terms of the fundamental string charge at the tip 
of the brane, which in the field theory corresponds to the density of massive fermions. In 
this way we can make contact with the "old fashioned" phase diagram of the GN model 
P^| . In this subsection we show how the corresponding effective potential can be obtained 
in the gravitational description. 

The D6 (or D6) probe brane action in the D4-background leads (after integration over 
the 5 4 ) to the 3-dimensional action 



5 6 = a fd 3 £ Ue-tyJ- det( 5afe + F ab ) ± k JA A F , (8.1) 
where a = TeR ^ 4 and k = The metric of the -D4 background at zero temperature 



is obtained from ( p.l|) with f(U) = 1. For a straight brane, representing the chirally 
symmetric phase, it is convenient to choose the coordinates £° = x ,^ 1 = x 1 and £ 2 = U, 
whereas for the curved brane we choose £ 2 = x 4 . It is not hard to show that on the straight 
brane the only regular solutions with finite action are pure gauge configurations. 
In the broken phase the probe brane action reads 



S 6 = a J d 3 £ uJ U' 2 + (J^J - 4A' + A'„ - kjd 3 £ (A_A' + - A + A'__) , (8.2) 

where ' denotes differentiation with respect to x 4 . The conserved (x 4 independent) quan- 
tities are 

P± = J]Y A ± =F kA± , 
U 4 

Pa = , = — ■ (8.3) 

Juv + (<l) 3 -4A> + A>_ 

In the absence of sources the only solutions with finite action are again the pure gauge 
configurations discussed in the last section. We will now consider the curved brane with 



- 22 - 



attached F-strings stretching along x 4 = from the brane tip at U = U$ to the singularity 
at U = 0. We will also assume that the F-string is smeared along x 1 . In field theory this 
corresponds to switching on a non-vanishing expectation value for the massive fermion in 
the broken phase. 

We construct this solution by taking the two different branches of source-free solutions 
from above, and matching them across the location of the source. Then, for a brane 
interpolating between x 4 = ±sl, we take 



A, 



3^P{-^9(U)) , x 4 <0 
, x 4 > 



f , x 4 < , 



iex P (-^(C/)) , x 4 >0 



y 3 



where 



and Uq is expressed through L by 

L f°° dy 



g(U) = / dy / , (8.5) 



(8.6) 



2 



The free parameters j and P4 are in turn fixed by matching this solution to the fundamental 
string source. To this end we add to the brane action 

S F = f dx°dUV- det h - f dx°A . (8.7) 

27TQ' J 27TQ' J 



Requiring the action to be stationary then yields 



2V 1 kj 



2ira' 

and 



Pr - m - , (8.9) 



where p = t-^ttt - Note that the above solution has no vector [/(l)-charge. Indeed using 



Q7.8D and (7.11) we we find fi = j and p = 0. In field theory language this means that 
the charge of the massive fermions has been compensated by a negative winding charge, 
kj, of the Goldstone boson eft. From our discussion in section 7.2. it is then clear that the 
asymptotic charge of this configuration can be adjusted to any value by a suitable axial 
gauge transformation. 7 



7 This gauge choice is actually the right one to compare the energy of this configuration with the straight 
branes with no gauge field on them. 
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For p < p max ~ 0.26 (^j^j the order parameter Uq (which determines the fermion mass 
through 0-7| )) is then a solution of the transcendental equation 

\m - 7p S{pm • <8 ' 10) 

where 

W^f vd-^-^ (8 ' n) 

For generic density there are two solutions for the brane profile (see Fig . These should 
correspond to the extrema of the Landau potential for Uq. To construct this potential 
we take the brane profile off-shell by giving up the jump condition (|8.9| ) but keeping the 
condition Q8.8| ) for the discontinuity of the gauge fields. We can think of this as displacing 
the tip of the brane "by hand" . To this end we write 

P 4 2 = R 3 UH (1 - x 2 ) , (8.12) 

where x S [0, 1] is our off-shell parameter. In this way all equations of motions will be 
satisfied apart from the jumping condition for U' . On-shell, from the .D6-brane point of 
view, the strength of the discontinuity in U' is a function of the strength of the external 
force, in our case, the fundamental string attached to it. 8 Uq is then the "order parameter" 
whose expectation value at fixed p (expressed in units of is obtained by extremizing 
the effective potential V(Uq). It turns out to be easier to determine the potential for x 
which is related to Uq by integrating the second conservation law in (IcO 



Ul' 2 {x) = ^— y/l 



2R 3 / 2 r s f dz 



X 



L J ^/z 8 — (1 — x 2 )z 3 

= ^ (i - x2) ~ 1/1 ° Bi -*44 ) - (8 - i3) 

Extremizing V(x) with respect to x should then imply the jumping condition ( |8.9| ), i.e. 
x = p/Uo(x). 

To continue, we consider the difference in energy between a curved brane with N 
smeared out fundamental strings and brane profile parameterized by x, and two straight 
branes with no gauge fields on them. Substitution of the gauge potentials and brane profile 
into (pJ] ), fl3,7| ) and (|8.2|) leads to the following expression for the energy difference: 

Ai [/ '" r) - pUq(x) + Uq(x) 2 I dz I - ~ ox = -z\ - U {xf I dzz 



2aVx 



P Uq{x) + U 2 {x){ 



z 5 




- (1 - x 2 )z 3 




1-x 2 




J dss~ ^(1 








-1/2 



.14) 



The first term on the r.h.s. of ( S.14| ) is the tension of the fundamental string (|3.2j). The 
remaining integrals in ( 8.14 ) are regularized incomplete beta functions. 

8 This is the analog of the external magnetic field for ferromagnetism. 
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The resulting potential is depicted in fig- 
ure [9|. The values of x for which F(p, x) at- 
tains its local minimum and maximum re- 
spectively correspond to the two solutions for 
Uq of ( |8.10| ). We will denote the former by 
x m i n and the latter by x max . We then see 
that the "second" solution (with larger value 
of x) has always higher energy and is there- 
fore unstable. 

For p < p c ~ 0.2 (in units of j^) the 
homogenous curved brane with fundamental 
strings attached to it has lower energy than 

the flat brane. Translating this back into field theory units, this gives 




Figure 8: Density, p as a function of x. 



10 



- 3 AiVA 
2vr 



5.15) 



For p c < p < Pmax — 0.26 the curved brane has higher energy than the flat brane (with no 
density on it). For p > p ma x there is no solution for the curved brane. The vanishing of 
the Landau potential at x = 1 may come as a surprise since for non-vanishing density the 
graphical representation (figure ||) of ( [3.10 ) does not give a solution with x = 1. 
However, inspection of the gauge field AF[P; x] 



configuration ( |8.4| ) shows that since P4 = 
for x = 1, the solutions for the gauge fields 
on the curved brane have vanishing support 
so that this configuration is in fact indistin- 
guishable from that of a curved brane with 
vanishing fundamental string charge and Uo, 
which is in turn a solution of (|8.10| ) . In other 
words pulling the tip of charged brane all the 
way down to the horizon one recovers the un- 
charged straight brane. 

The picture that arises here appears to 
be in qualitative agreement with the mean 




Figure 9: Effective potential for x for p = 0.1, 
p = p c = 0.2 and p — p max — .26 (in units of 
R 3 /L 2 ). 



field approximation (see e.g. [31]) assuming one assigns a unit charge to the fundamental 
string. Then our effective potential for the mean field Uo predicts a phase transition to the 
chirally symmetric phase at p = p c in agreement with the (old fashioned) phase diagram 
for the GN model. 

Finally one can obtain the mass M of the "mean field fermion" as 



2vra'M = 



dE, 
dp 



\p=0 ) 



T U + — \x=0 \ 
5 b ox 



dss~ 



(1-s) 



-1/2 
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(8.16) 



This is not the physical fermion mass in general. Rather it is the mass of a delocalized 
fermion obtained by assuming that single particle states are momentum eigenstates (which 
would be correct for a free theory). If these massive fermions were non- interacting we 
should get 2WM = U . 



8.2 Stringy NJL model 

We now turn to the stringy NJL model. Again we will consider the model at finite density 



but zero temperature for simplicity. In this case the matching conditions (3.12) for a 
smeared fundamental string source become 9 



Q = 7^7 = P^R* , 



2ira' 



and 



P 



U \ 



Q "77 



R J 



5.17) 



(8.18) 



Upon substitution of fl8.17|) the later leads to (|3.13|) at T = 0. We could now define the off- 
shell Landau potential by relaxing the constraint (|3.13j ) for P4 while keeping ( 8.17 ) as we 



did in the case for the D6-brane. However, unlike for the D6-brane where Q did not appear 
in the DBI action and the expression for U, this ansatz leads to complicated expressions in 
which Uq is only implicitly determined, so that one has to solve the corresponding system 
numerically as in section 4. For the sake of computational simplicity we will therefore follow 
a different strategy, in which both ( |3.17 ) and ( |8,18j ) are taken off-shell simultaneously by 
introducing a single parameter x. Since both ( 8.17| ) and (3.18) are equations of motions for 



the combined system of probe brane and fundamental strings, this should be a consistent 
off-shell definition, at least near the critical points. Concretely we set 



x 



Q 2 



K 2 R 3 Ul(a 



and 



P A {x) = kU$(x) 



1.19) 



where Uq(x) is determined by integrating (3.13), i.e. 



U (x) 



4R S 



with 



{9{x)f , 



dz 



11 



Z A + z°x 



6^2 



1.20) 



121) 



Note that if (8.17) is not satisfied then P4 as defined in (8.18) is taken off-shell as well, since 



the charge no longer drops out in (8.18j). For a generic value of x the jumping conditions are 



9 We absorb (3 into the definition of the charge. 
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satisfied for some charge density which is different from that sourced by the fundamental 
string charge p. To see this we plot 



^ = , (8-22) 

against x (see figure [n]). For p < p max — (0.47) 5 (i?2 /L) 5 there are two solutions for 
fixed fundamental string charge p. To see which solution has lower energy we compute the 
Landau free energy off-shell as a function of x at fixed p. The actual solution is obtained 
by minimizing the Landau free energy with respect to x. In order to get the free energy 
we take the Legendre transform of the brane and string action, as explained in section 4, 

oo 

F(Q) = 2 J L DBI (A, A') + S F + 2QA(oo) . (8.23) 

Uo 

The factor of 2 takes into account both legs of the curved brane. For a given charge Q(x) 
on the curved brane we have 

OO 3 

A curved (x,oo) =xU (x) f ; dZZ2 , (8.24) 

l 

while for a straight brane with charge p we have 

A straight^ Qo) = p \ ( g 25) 



7T 

The free energy difference between the curved brane with charge Q{x) and the straight 
brane with charge p is then given by 10 

AF(p,x) = -S F (p,U (x)) + 2 P kR^(A c (x,oo) - A s (p,oo)) 

- 2Sh BI {A'(x)) + 2Sh BI (A'(p)) . (8.26) 

Putting it all together we get 

— = pUo(x) +pxUo{x) I /a = ==-p5. 

2/ti?2 



.27) 




where the order of the different terms in ( |8.27 ) is the same as in ( ^.26 ). 



The remaining integrals are easily evaluated numerically. Figure 11 showing the ef- 



fective potential for x reveals a picture that is in quantitative agreement with figure 10 



10 Since the Legendre transform of the grand canonical potential is well defined only on-shell (i.e. when 
Q is the physical charge of the system) there is some ambiguity in defining F(p, x). Here we take the point 



of view where the Q entering in (S.23) is the physical charge. 
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Figure 10: Density p as function of x as in (8.22) 



and ( |8.22| ) in particular. For given charge density p with p < p ma x the effective potential 
has a local minimum, which corresponds to the smaller value of the 2-possible solutions of 



.22), and a local maximum which corresponds to the larger value of x. The actual value 



of 

Pmax obtained by analyzing the effective potential is identical with ( |8.22| ) (see figure 
11). Another feature of the effective potential is that the energy of the curved brane with 
finite charge density exceeds that of the straight branes for p > p c with p c ~ (0.23) 5 . The 
Landau potential predicts a first order transition to the chirally symmetric phase. 



As explained in section 5 the mean field 
approximation discussed here does not yield 
the correct phase diagram at large N c . The 
physically realized phase at non- vanishing den- 
sity is rather a inhomogeneous mixture of 
phases of uncharged curved branes and charged 
straight branes. Although we will not repeat 
the analysis here, it is not hard to show that 
the free energy of the mixed phase is lower 
than the minimum of the Landau free en- 
ergy of the pure phase. The corresponding 
phase diagram is identical with that obtained 
in section 5 for zero temperature. 



AF [p, x] 



-5-1 




Figure 11: A_F(p, a;) as a function of x for p = 
(0.2) 5 , p = Pc = (0.23) 5 and p = p max = (0.24) 5 
in units of (2 J R 3 / 2 /L) 5 . 



We conclude this section by noting that the Landau potential does not vanish for 
x —* oo unless p = 0. The fall-off of AF for x —* oo is a consequence of the unphysical 
relaxation of the matching condition for the charge in ( |B.17 ). The corresponding Landau 
potential is unphysical for large x. In particular, the fall-off at large x should not be 
interpreted as an instability of the configuration corresponding to the local minimum. 
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9. Discussion 

To conclude, we briefly summarize the main results of this paper. We have mapped out the 
finite density/temperature phase structures of the stringy NJL and GN models at strong 
coupling by studying probe branes in the near horizon geometry of D4-branes. In the 
stringy NJL model an important role was played by solutions with fundamental strings 
having one end on the probe branes and the other end disappearing through the D4-brane 
horizon. These fundamental strings carry baryon number in the phase with broken chiral 
symmetry. We also reviewed the corresponding phase diagrams in the local field theory 
cousins of these models, and compared with the strong coupling supergravity results. For 
the stringy NJL model we saw good qualitative agreement in the phase structure at low 
temperatures, except at the very highest densities. On the other hand, at high temperatures 
there is a qualitative difference between the two sides, the chiral phase transition being 
respectively first and second order in the supergravity and field theory models. For the 
stringy GN model there is a precise agreement between the two sides, but in a somewhat 
trivial way. In both supergravity and field theory, the charge density is simply carried by a 
free boson. In supergravity the free boson corresponds to pure gauge modes on the probe 
brane, and in field theory it appears via bosonization. 

In our study of the phase structure we encountered various stable, unstable, and 
metastable phases. For instance, when we looked for curved brane solutions with attached 
fundamental strings, we found that the equations of motion admitted multiple brane pro- 
files with a given charge density. A convenient way to keep track of these different solutions 
is via an off-shell Landau potential. By taking the probe branes off-shell in a particular way, 
we can interpolate between the various solutions, and the absolute minimum of the Landau 
potential determines the stable phase. We illustrated this procedure in the simplified case 
of zero temperature. 

As discussed in the introduction, the long term motivation for the analysis considered 
here is to eventually apply stringy methods to learn about finite density QCD. Since NJL 
type models are one of the main phenomenological tools in this subject, studying their 
string theory incarnations is potentially useful. The main challenge ahead is to find ways 
to modify the stringy models so as to make them more closely resemble physical QCD. 
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